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Q ' Abstract 
(N ■ 

The Gaudin models based on the face-type ehiptic quantum groups and the 
[ XYZ Gaudin models are studied. The Gaudin model Hamiltonians are con- 

structed and are diagonalized by using the algebraic Bethe ansatz method. The 
■ corresponding face-type Knizhnik-Zamolodchikov equations and their solutions 

^ ^ are given. 
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I Introduction 

In [|I|], Gaudin proposed a quantum integrable model describing spin 1/2 particles 
with long-range interactions. The Gaudin type models played an important role in 
^ i establishing the integrability of the Seiberg-Witten theory [|, and were used as a 
QQ I testing ground for ideas such as the functional Bethe ansatz and general procedure of 
' separation of variables 0, ^, |^. They also have direct applications in condensed matter 

physics. 

The algebra associated with the model proposed by Gaudin is the Lie algebra 5/2- 
Gaudin's work was later generalized by several authors J^, §|, |§] and integrable Gaudin 
I models based on other Lie algebras were constructed. The XYZ Gaudin model was 
counstructed and solved in |^ by means of the algebraic Bethe ansatz. The boundary 
^ ■ Gaudin magnet associated with the spin 1/2 representation of SI2 was investigated by 

> : Hikami |0| using Sklyanin's boundary quantum inverse scattering method. 
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^ i The Knizhnik-Zamolodchikov (KZ) equations were first proposed as a set of differ- 

^ I ential equations satisfied by correlation functions of the Wess-Zumino-Witten models 
|n]. The relation between the Gaudin magnets and the KZ equations has been studied 



in many papers |T^, 0, |1J, |T5|. In |TJ] and [|T^], Hikami gave an integral representation 
for solutions of the KZ equations by using the results of the periodic and boundary 
Gaudin models. 

In the present paper, we construct elliptic Gaudin models based on the face-type 



elliptic quantum group Er^riish) [|T6[ and boundary elliptic quantum group BEr^nish) 
f[% and the boundary XYZ Gaudin models. We diagonalize them by means of the 
algebraic Bethe ansatz method. Moreover we construct the face-type elliptic KZ equa- 
tions and their solutions using the off-shell Bethe ansatz equations of the face-type 
elliptic Gaudin models. 

This paper is organized as follows. In section 2, we review the elliptic quantum 
group Et-^^{s12) and the boundary elliptic quantum group BET-r,{sl2). Then in section 
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3, we construct the corresponding elliptic quantum Gaudin models. Solutions of the 
KZ equations based on these Gaudin magnets are also given. In section 4, we give the 
Hamiltonians, eigenvalues and Bethe ansatz equations of the boundary XYZ Gaudin 
model. In the last section, we present some discussions. 

II Elliptic quantum groups and Bethe ansatz 

In this section we review the elliptic quantum group i?T-_^(s/2) [0, and the boundary 



elliptic quantum group BEt-.^{s12) [|TT1- These elliptic quantum groups are algebraic 
structures underlying the (boundary) dynamical Yang-Baxter equation in statistical 
mechanics and the KZ equations on a torus. 

II. 1 £',-^(5/2) and its Bethe ansatz 

Let h be the generator of an one-dimensional commutative Lie algebra, and a{X,w), 
b{X,w), c{X,w), d{X,w) be the elements of the matrix L{X,w). Then Er^rjish) is 
generated by meromorphic functions of the variable h and the elements of L(A, w) with 
non- commutative entries, subject to the dynamical Yang-Baxter equation (YBE) 



_ 27]h,wi2)L^^\X,wi)L^^\X - 2r]h^^\w2) 
= L^^\X,W2)L^^\X-2T^h^^\w,)R^^^\X,w,2). (H.l) 

Here Wij = Wi ~ Wj with Wi being spectral parameters, 77 is the crossing parameter; 
L*^')(A — 2rih^^\ Wi) acts on the z-th auxiliary space and h^^^ on the j -th auxiliary space 
^ = by a diagonal matrix diag{l, —1). R^^'^\X~2rih, W12) acts on the space C^(S>C^ 
with h acting on the quantum space. The i?-matrix is given by 

R{X,w) = Eo,o®^o,o + ^i,i®^i,i + tt(A,w)Eo,o®^i,i + «(-A,w)Ei,i®Eo,o 

+ P{X, w)Eo,i ® ^1,0 + Pi-X, w)E^fl ® Eo,i (II.2) 

with 

e{w)6{X + 2^) e{X + w)e{2r^) 

"^^'^) = ^(«;-2r,MA)' = ^(«; - 2r,MA) ' 

d{X) = 9{X,t) = - ^ g7ri(j+l/2)2r+27ri(j+l/2)(A+l/2)_ 

The commutation relations between the generators of Er^r^ish) are defined by 
eq. (|II.l| ), and the commutation relations of the generators with some functions /(A, h), 
giX,h) are PI 



/(A, h)f{X, h) = g{X, h)f{X, h), /(A - 2r], h)a{X, w) = a(A, ^)/(A, h), 
/(A + 2r7, h)d{X, w) = rf(A, w)f{X, h), /(A + 2r],h + 2)6(A, w) = 6(A, w)f{X, h), 
/(A ~2r],h~ 2)c(A, w) = c(A, w)f{X, h). (ILS) 
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For an even number A > and a complex number z, we can define an evaluation 
module Va{z) of £^,-^^(5/2)- Let e^, k G Z>o be a set of bases of V\{z). The action is 
defined by 

fih)ek = /(A - 2k)ek, a{\, w)ek = g{a, X, w, k)ek, 

b{X, w)ek = g{b, A, w, k)ek+i, d{X, w)ek = g{d, X, w, k)ek, 

c{X,w)ek = g{c,X,w,k)ek-i, (II.4) 



with 

g{a,X,w,k) 
g{b,X,w,k) 
g{c,X,w,k) 



9(z-w + (A + l- 2k)r])9{X + 2kr]) 
e{z - w + {A + l)7])e{X) ' 
^(-A + z-w + {A-l- 2k)r])9{27]) 
e{z ~ w + {A + l)7])e{X) ' 
^(-A -z + w + {A + l- 2k)ri)9{2{A + 1 - k)ri)9{2kri) 
e{z - w + {A + l)r])e{X)e{2r]) ' 



e{z-w+{-A + l + 2k)r^)0{X - 2(A - k)^) 

^^^'^'"'^^ = e{z-w + iA + i),)9ix) • ^"-'^ 

For any finite-dimensional module W = V\ of Et-^^{s12), the transfer matrix can be 
defined by 

t(A, w) = troL{X, w)) = a(A, w) + d{X, w). (II.6) 

Then t{X,w) preserves the space H = Fun{W)[0] of functions with values in the zero 
weight space and commutes pairwise on H, i.e. t{X,w)t{X,u) = t{X,u)t{X,w) on H. 

Let W = Vai(-2i) ^ ■ ■ ■ ® Vaj^{zn) be a tensor product of evaluation modules of 
Er^riish) and let A = Ai + ■ ■ ■ + A^. Then, the highest weight state of the space H 
with W[A] = Cvq obeys the following highest weight condition, 

c{X,w)vq = 0, a{X,w)vQ = vq, 

^ , 9{X- 2r]A) 9{w - z. - (-A + l)r/) 

Under the framework of the algebraic Bethe ansatz, the Bethe state is defined by 

■ ■ ■ , tAf) = b{ti)b{t2) ■ ■ ■ b{tM)v = b{t,)b{t2) ■ ■ ■ b{tM)g{X)Vo, (IL8) 

where g{X) 7^ is a meromorphic function. 

According to Felder and Varchenko |]16|, W is a. highest weight module of 
Er,r,{sl2) with A = 2M G 2Z>o. Let v{X) = n|li ^(A - 2r/j). Then applying the 
transfer matrix ( |II.6|) to the Bethe state gives rise to 

t(A,w)$(ti,---,tAf) 
= {a{X,w) + d{X,w))<l>{ti,---,tM) 

M 

= e{X,w)^{ti,- ■ ■ ,tM) + ^xFa^aitir ■ ■ ,'ta-l,W,ta+l,tM), (11.9) 

j=a 
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where Fa is a function of ta, x = d{X + ta — w)6{2r])/{6{ta — w)6{X — 2r])) and 

^a(ti, ■ ■ ■ , ta-i, w, ta+1, Im) = b{X, h) ■ ■ ■ 6(A, ta-i)b{X, w)b{X, t^+i) ■ ■ ■ 6(A, tM)v. 

It follows that the Bethe state is an eigenstate of the transfer matrix t{X, w) if = 0. 
Therefore the eigenvalue of the transfer matrix is 

' Ik ^(^« - iii ^(^- - M ^(«^ - - (A^c + 

(11.10) 

with ta determined by the Bethe ansatz equations, 

p _ 1 n ^^^^ ~ ~ TT ~ ~ + ^^"^^ - n mill 

^}a(^{tp-ta + 2,^)l}^ e{ta-Z,-{K + l)v) ' ^' ' 

II. 2 BEr^rjish) and its Bethe ansatz 

The boundary elliptic quantum group BEr^r^ish) is generated by h and the elements 
of the matrix C{X,w) G End(C^), A{X,w), B{X,w), C{X,w) and V{X,w) with non- 
commutative entries, subject to the relations 

i?2i(A, wi - W2)Ci{X - 2rih^^\wi)Ri2{X, wi + ws) A(A - 2rih^^\w2) 
= ^(A - 2r7/iW, w;2)i?2i(A, wi + w;2)/:i(A - 277/^(2), w;i)i?i2(A, - ^2). 

(11.12) 

Let V\{z) be an evaluation module of BEt.,^{s12) with the basis e^, k G Z>o . This 
module can be identified with the original representation space of Et-^^{s12)- On the 
module Va(-2), the representation of ^£'^^^(5/2) is given by 

f{h)eu = f{A-2k)ek, (11.13) 

A{X,w)ek = {Y{X + 2ri,—w)g{a,X + 2ri,w,k)g{d,X + 2ri,—w,k — l) 

- K{X - 2r]{A -2k + 2),w)22Y{X + 2r], -w) 

X g{c, X + 2?7, w, k)g{a, X + 2r], —w, k — l)}ek, (11-14) 

B{X,w)ek = {Y{X + 2r],-w)g{b, X + 2r],w,k)g{d, X + 27], -w,k) 

- K{X - 27]{A - 2A;), w;)22l'(A + 2r], -w) 

X g{d, X + 2ri, w, k)g{b, X + 2ri, ~w, k)}ek+i, (11.15) 

C{X,w)ek = {K{X-2ri{A-2k + 2),w)22YiX-2r],-w)X{X-2ri,k- 1) 
X g{c, X — 2?7, w, k)g{a, X — 21], —w, k) 

- y(A - 2r], -w)X{X -2r],k- 1) 

X g{a, X - 2t], w, k)g{c, X - 27], -w, k)}ek-i, (11.16) 
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V{X,w)ek = {K{X-2r]{A-2k),w)22Y{X-2r],-w)X{X~2ri,k) 
X g{d, X — 2?7, w, k)g{a, X — 2ri, —w, k + 1) 
-Y{X-27],-w)X{X-2r],k) 

X g{b, X - 2r], w, k)g{c, X-2r],w,k + l)}efc, (11.17) 



where ^ is an arbitrary parameter, g{i, X,w,k), i = a, b, c, d, are given by (|lT 

e{X-2{A + l-2k)7]) 



X{X,k) 
YiX,w) 



e{X-2{A-l-2k)r])' 

e{z - w + {A + l)r])e{X) 



e{z-w-{A + l)v)0{X - 2(A + 1)7]) ' 
and K{X,w)22 is the (2,2) element of the K-matrix 

which satisfies the boundary dynamical YBE: 

i?2l(A, Wi - W2)K,{X - 27]h^^\wi)Ri2{X, Wi + W2)K2{X - 27]h^^\w2) 

= K2{X - 2'qh^^\w2)R2i{\wi + W2)Ki{X - 2rih'^^\w{)Ri2{X, Wi - W2). 

(11.19) 

Similar to the Er^n^sh) case, for any module = Va of BEr,t^{sl2), the boundary 
transfer matrix is defined by 

t\X, w) = tToK+{X, w)C{X, w), (11.20) 

where 

is a diagonal solution to the dual boundary dynamical YBE, 

R2i{X, -wi + W2)Kt{X - 27]h^^\wi)Ri2{X, -wi - W2 + 4:r])K+{X - 2r]h^^\ W2) 
= K+{X - 2r]h^^\ W2)R2i{\ -wi -W2 + ^v)Kt{\ - 2r]h^^\ Wi) 

X Ri2iX, -wi + W2) (11.22) 

with 

R{X,w) = ^0,0 ® ^0,0 + ^1,1® ^1,1 + a(A,w)Eo,oO ^1,1 + a(-A,w)^i,i 0^0,0 

(11.23) 

The transfer matrix ( [II. 201) preserves the space H = Fun{W)[0] of functions with 
values in the zero weight space W[0]. Let W = Vai(-2i) ® ■ ■ -Vaj^^zn) be a tensor 
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V(X,w)vo = 



product of evaluation modules of BEr,n{sl2), and let A = Ai + ■ ■ ■ + A^. The highest 
weight vector vq of H with 1^[A] = Cvq satisfies the highest weight condition: 

C{X,w)vo = 0, A{X,w)vo = vo, 

e{2w)e{w - ^ + A - 2(A + l)r])e{w + ^-2r]) 
e(2w - 2r])e{w + e - (A - 2A7]))e{w - 
SL e{z, +W + (A, - l)v)9iz, -w- (A, - 

where T>{X,w) is determined from the relation: V{X,w) = 9{X)/6{X — 2ri)T>{X,w) + 
P{X — 2r], 2w)A{X,w). The space H = Fun(W) is spanned by the Bethe state 

^'(^1, • • • , ^m) = B{h) ■ ■ ■ B{tMW (11.25) 
Applying the boundary transfer matrix to the Bethe state, one obtains 

M 
a=l 

(11.26) 

where 



^(277)^(2^,) 

X 



e{X + 2rj)9{2t^-2rj) 



0{2w - Ar])e{w + ta - X) e{w - ta - X - 2r)) 
9{2w -2rj)9{w + ta-2rj) 9{w - t^) 



. . _ yr 9{w + ta)9{w - + 2r]) 9{w -2ri- i)9{w - A - 2r/ + Q 
' iii 0{w + t„ - 2r^)9{w - t J ^ 9{^v + A - + ^ 

^(2w)^(u; - 2r/ + + A - 2(A + l)r/ - Q 

^ 9{2w - Ari)9{w - ^)9{w - A + 2A77 + 

^ 9{w + tg- 2r])9{w - tg) 

^ n ^(^^ + ^ + (A^-lhM^^-^-(A/-l)r]) , . 

l}^9{zi+w-{Ai + l)7^)9{zi-w + {Ai + l)r]y ^ ' > 

and obeys the Bethe ansatz equations: 

^i^^^ 9{to, + i;3)^(ia - + 2ri) 

^ 9{t^-2ri-09{t^-2ri-X + 

0{ta + X-O0{ta+O 

9{ta - 277 + i)9{t^ - 2(A + 1)77 + A - e) 

^ e(ia - A + 2A77 + 0^(^a - 

^(^; + + (A/ - l)r/)^(^^ - - (A; - 1)77) 
tk ^i^i + (A; + l)v)0{zi -ta + {Ai + 1)77) 
= 0. (11.28) 
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III Gaudin magnets and KZ equations based on the 
elliptic quantum groups 

III.l Periodic Gaudin magnet 

The Gaudin model was proposed as a quantum integrable system with long-range 
interactions [0]. Such a system can be solved by using the algebraic Bethe ansatz 
method. By taking the quasi-classical limit rj of the transfer matrix of the six- 
vertex model, Hikami gave the Hamiltonian of the XX Z Gaudin model |]T3| . 

Foe our case, we expand the transfer matrix ( [1I.2CI| ) around the point r] = with 
the parameter w = Zj, to obtain 



t(A, w = z,) = + VH, + 0{r,')). 



(III.l) 



One can prove [ifj,iffc]=0 since the first term on the right hand side of ( [1II.1| ) is a 
constant. Therefore, Hj gives the Hamiltonian associated with Er^nish)- From ( [II. 6|) , 
we find 



dt{\, Zj 



di] 



ri=0 



N 



J2 [Wiizi, z,)EfE+ + W2{zi, z,)EtEj - 2h,dx 

1 1 

+ -ip{X){Aihj + hi) + -(p{zj - zi){hihj - Ai) 



where ip{x) = 9'{x)/9{x) 



Wi{x,y) 
W2{x,y) 



9{-\ + x-y)e'{0) 

o{x-y)e{X) ■ 
e{-x-x + y)e'{<d) 
e{x-y)e{\) ■ 



and h, E^ are defined by 



hiCk = (A; - 2k)ek, 

EfCk = Cfc+i, E^Ck = k{Ai + 1 - k)ek- 
We can check that h, E"^ satisfy the commutation relations 

[K, Ef] = ±25,,Ef, [Et, Er] = 6.,h,. 

The eigenvalues ej of Hj can be extracted from 

2 



e(A, w = z,) = ^^^(1 + V^j + Oiv')), 



(111.2) 



(III.3) 



(111.4) 



(III.5) 



(III.6) 
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M N 

= -2Aj vita - Zj) + (1 - Aj) - zi). 



(III.7) 



a=l 



And the constraints for the eigenvalues are fa = 0, a = 1,2, ■ ■ ■ , M, with 

M N 

fa = ^4¥p(t/3-t«) + 2^A,¥?(t„-2;i). 

(3^a 1=1 



(111.8) 



We now define the following equivalent Hamiltonian Tij of the Gaudin model by 
shifting Hj by some constants, 



N 



n, = H, - (1 - A,) Y ^iVizj - zi) 



(111.9) 



The corresponding eigenvalues are 



M 



(111.10) 



Then, from ( |11.9|) , we obtain the so-called off-shell Bethe ansatz equations: 

M 

UjCj) = £j(f) + J2 fcWl{Zj, ta)E-(j)a, 
a=l 

where 0, (pa are Bethe states for the Gaudin model given by 

10 >, 



(111.11) 



^ db{X,ta) 



ri=0 



M 



- n 



db{X,t 



dr] 



10 > . 



(111.12) 



??=o 



By using the representation ( [11. 4|) of Et^^^sI^), we derive the precise form of the Bethe 
states 

M N 

^ = l[J2W^i^>^^ta)E^\0>, 

a=l k=l 
M N 

^j^a k=l 



III.2 Elliptic KZ equation 

As integrable differential equations, the KZ equations take the form 

Vj* = for j = l,2,---,Ar, 



(111.13) 



(111.14) 
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where the differential operators Vj are defined by the Gaudin model Hamiltonian Ti.j 

d 



(III.15) 

with K being an arbitrary parameter. Substituting (|1II.2D into (|III.15|) , we can check 

[V,,V,.] =0, (III.16) 

which ensures the integrability of the KZ equations. 

The function \E'(-2) can be constructed by the hypergeometric function t) which 
is a solution of the following equations 

^ f 

Otoe 

The hypergeometric function xi^z^ t) can be written as 

M N 



(III.17) 



x{z,t) = n ma - 1,)]-'/^ n um^ - ^«)]'^^/^- 

/3<a a=l j=l 



(III.18) 



With the help of x{z,t), we obtain 



C JC 



dti ■ ■ ■ dtMXi't, z)(l){t, z) 



(III.19) 



where the integration path C is a closed contour in the Riemann surface such that the 
integrand resumes its initial value. 

Substituting the expressions of Vj and "^{z) into (|1II.14| ), we can show that the KZ 
equations are satisfied. The proof is as follows 



9 T, > 



/ ■ ■ ■ / dti--- dtM \i^^<P + I 

JC JC \ OZj OZj J 

■ ■ ■ rfti ■ ■ ■ dtM + ^^^^ 

f ■ ■ ■ f dti - ■ ■ dt 

JC JC 
M 



d 



a=l Ol^a 



f ■ ■ ■ f dti - ■ ■ dtM 

JC JC 



ic JC 

n.^iz), 



'HjX<P - {xWi{zj,to,)E]' 

a=l OI^O' 



(III.20) 
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III. 3 Boundary Gaudin magnet and KZ equation 

The Gaudin model Hamiltonian for BE^-rjlsh) is found to be the form 



H'{X,z^) = 6{X,zj,0 



N 



k=l 



n 

+ E [W{{z,,z^)E,El + W^{z,,z,)EtEj 
+ ip{zk - Zj){hkhj - Afc) - ip{\){hk - Akhj)] , 



(III.21) 



where E^, h and (p are the same as in the periodic case, K22 and K22 are (2, 2)-elements 
of K and respectively, and 



S{X,Zj,^) 



e{-x + x±y)e'io) 
eix + x±y)e'iQ) 

9{x±y)9iX) ' 
(1 + 2hj)ih^ - l)[^iz, - A + - V^(^, + A - 0] 
+ {h, - 1)[2^(A) + ^{z, - - V^(% + 0] 



+ h,{A,-l) 



(III.22) 



By ( [11.271 ) we find the eigenvalues of the boundary Gaudin model Hamiltonian Hj, 

de{X, w = Zj) 



dr] 



r)=0 



M 



E 2Aj[v5(zj + ta) + ^{zj - ta)] - (1 - A^j)ip{2zj) 



a=l 



N 



+ E - ^3)Hzi + Zj) - ifiizi - Zj)] 

- (1 - A,)b(% + + ¥^{^j - + <^(% -X + + <^(% + A - 0] 

(III.23) 

with the constraints = 0, a = 1, 2, ■ ■ ■ , M, where 

= 2(^(t,-0 + 2<^(t„ + + 2(l + A)[<^(ta-A + + <^(ta + A-0] 

A/ Af 
+ 4 E Hta-tf3)+ipit^+tp)]-J2'^Mvizi+t^)-ipizi-t^)]. 
P=l,^a 1=1 

(III.24) 
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The Bethe states for the boundary Gaudin model are then given by 



0^ 



M 

n 



N 



a=l 

M 

n 



1=1 



N 



.1=1 



(III.25) 



As in the previous subsection, we define the following equivalent Hamiltonian Hj 
of the boundary Gaudin model: 



AT 



- (1 - Aj)[zj + + <p{zj - + <P{zj - A + + <P{zj + A - 0]- 



The corresponding eigenvalues are given by 



M 



E] = (1 - A|)^(2z,) - 2A,[^(^, + to) + ^{z, - 



(III.26) 



(III.27) 



a=l 



and the off-shell Bethe ansatz equations read 

M 



n)<t>' = £:>^ +Y.fi {Wii^u ta) - W+izi, t^)) E-<f>l (III.28) 



a=l 



Thus the KZ equations based on the boundary Gaudin model are given by 

V.*'' = for i-l,2,--.,Ar, 



where the differential operators Vj are defined by 



(III.29) 



(111.30) 



with K being an arbitrary parameter. To ensure the integrability of the KZ equation, 
we impose 

[V,-,Vfe]=0, (III.31) 

This requires ^ oo. With this condition, the Hamiltonian and the Bethe states of 
the boundary Gaudin model become 

N 



= -^{\z,,i)-Y.[wtE^Ej + W^EtE. 



k=l 



- (p{zk + Zj){hkhj - Afc) + (p{X){2Ak - hk - Akhj)] 

n 

- E [w,-E,El + W,-E^Er 

k=l,^3 
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+ (p{zk - Zj){hkhj - Ak) - (p{X){hk - Akhj)] , 



M 



a=l 11=1 



N 



M 

= n 



N 



J2{W,-{zi,tf,)-W,+ {zi,tf,))Ef 



(III.32) 



where 



N 



A(A) 



The function can be computed by using the following hypergeometric type in- 
tegral 

^\z) =£...£ dti •• • dtMx\t, z)(t>\t, z), (III.33) 

where the integration path C is a closed contour in the Riemann surface, similar to 
that described in the previous subsection. The hypergeometric kernel x^(^, z) is given 

by 



N 



N M 



j=l j=l a=l 

M 

X lli^{ta-h)9{ta + tfs)Y/\ 
a</3 

One can check that the kernel x''(t, z) satisfies the following equations 

^^x' = Six, 



(III.34) 



dzj 



Moreover we can show that satisfies the KZ equations, 

V^^\z) = 0. 

The proof is as follows. 

K^-^iz) = <f ■■■ (f dti---dtM(K^(l)''+Kx'^ 
OZj Jc Jc \ OZj OZj 

dti---dtM W"H](j)^ 



(III.35) 



(III.36) 



c JC 
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M 



M o 



(iti ■ ■ ■ dt][.f 



b 

3 "^a 



- E Kr(^.'^") - Ejci>i 



(III.37) 



IV XY Z Gaudin magnets 

The XYZ chain or eight-vertex model occupies an important place in the study of 
integrable systems pO[. The periodic XYZ Gaudin model has been constructed and 
diagonalized in |^ by means of the algebraic Bethe ansatz method. The corresponding 
KZ equations and their solutions were given in Here we write down a more 

explicit formula for the Bethe states for the periodic case. Moreover, we present the 
Hamiltonian, eigenvalues and (off-shell) Bethe ansatz equations for the XY Z Gaudin 
magnet with an integrable boundary. 

The Boltzmann weights of the eight-vertex model are given by 



R{w) 



Eqo ® Eqo + Ell 
snfw' 



+ 



Ell + K,sn{ri)sn{w){Eoi ® Eqi + Eio 

sn{ri) 



sn{w + 7]) 



{Eoo ® Ell + Ell ® 



00) 



sn{w + T]) 



Eio) 

{Eoi^Eio + Eio^Eoi), 
(IV.l) 



where sn{w) = sn{w, k) is the Jacobi elhptic function of modulus k, < n < 1. Define 
the corresponding transfer matrix as 

t(w) = troRoi{w - zi) ■ --RoNiw - zn) (IV.2) 

for the periodic case, and as 

t\w) = troK^iw,ORoiiw-zi)---RoNiw-ZN)K{w,^) 

xRon{w + zn)---Roi{w + zi) (IV.3) 

for the boundary case, where K{u,^) = diag(sn(^ + 'u),sn(^ — 'u))/sn(^) and 
K'^lu,^) = K{~u — ?7, ^) with ^ being a parameter. 

Then, by using the same method as in the previous sections, one can obtain the 
Hamiltonians of the XYZ Gaudin models. 



(i) The periodic case: 



N 

E 



[1 + Ksn^{zj - Zk)crja^ + (1 ^ Ksn^{zj - Zk)ajal 



+cn{zj - Zk)dn{zj - Zk){(Tj(Tl - 1] 



(IV.4) 
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where and throughout this paper, cn(n)dn(u) = dsn{u)/du. 
(ii)The open boundary case: 



H 



1 



2sn(^ - Zj) 
1 

" 2sn(e + Zj 

N 



cn(^ - Zj)dn{^ - zj) - 
cn(^ + Zj)dn{^ + Zj) + 



sn{2zj) 
sn(^ - Zj 
sn{2zj) 



K - 1) 



T — = — 

fcfl^j sn(2;j - Zk) 



+c\i{zj - Zk)dn{zj - Zk) 



1 



N 

^ sn(2;,- + Zk) 



KSr?{Zj + Zk) 



( sn(^ + Zj) 



\^\\.{i-Zj) 



0' ^+ 



sn(^ - Zj) 



sn(^ - z^ 



sn(^ + Zj 



(IV.5) 



The Hamiltonians ( [IV .41) and ( |IV.5| ) can be diagonahzed by using the method sim- 
ilar to the previous section. The results are as follows. 

(i) The periodic case. 

The eigenvalues of (|IV.4| ) are found to be 



M 



Ej = -2 ^ [VA^Zj - to) + ^i{zj - to)] 



(IV.6) 



a=l 



with the constraints fa = ^i = 2, ■ ■ ■ , M, where 



M N 
/a = -4 [Vi{tp-ta)+'^l{tl3-t^)]-2Y,[^i{to,- Zj)+'^l{t^- Zj)]. (IV.7) 

/3=l,^a j=l 

Here ^^{z) = Q'{z)/Q{z) and y:>i{z) = H'{z)/ H{z). Q{z) and H{z) are the Jacobi 
theta-functions and they satisfy the relations 

sn{z) = H{z)/{Vke{z)), sn\z) + cn\z) = 1, khn\z) + dn\z) = 1. 

We can write down the off-shell Bethe ansatz equations 

1 



AI 



Hj(j) = E,(j)+J2xf> 



where 



X 



e(t„ - Zj)H{ta - Zj 

2g'{0)g{^ + t^-z,-l) 
9{ta - zM'i^ - i) 



-iV(t)0, 



(IV.8) 
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with s,t being parameters and g{u) = H{u)Q{u 

N{t) = c{s,t) 



Qit)H{t) -H{t)H{t) 
0(t)0(t) -Q{t)H{t) 



The Bethe states (p and 0„ are given by 

M N 



M N 1 

n Etttt Trni -Mt,tp,z^)\Q>, (IV.9) 



where 



with 



Mil Mi2 

M21 M22 



Mil = 2aQ{ta - Zj)H{t + - Zj)[Q'{t + t„ - Zj)H{t + t„ - Zj) 

- 0(t + - + ta- Zj)] 

+ e(t + ta - Zj)H{t + ta- Zj)[e\ta - Zj)H{ta - Zj) 

-Q{t^-Zj)H'{to^-Zj)l 
H'(0) 

M12 = -Q^[0'(^ + ta- Zj)H\t^ - Zj) - e\t^ - Zj)H\t + t„ - Zj)l 

H'(0) 

M21 = -Q^[0'(^ + tc.- Zj)Q\t^ - Zj) - + - Zj)H''{t^ - Zj)\ 

M22 = 2ae(t« - Zj)H{t + t„ - Zj)[Q'{t + t^- Zj)H{t + t^- zj) 

- e(t + 1„ - Zj)H'{t + t^- Zj)] 

- e(t + 1« - 2j)i/(t + 1« - ^i)[e'(t« - Zj)H{t^ - Zj) 

-Q{t^- Zj)H'{t^- Zj)]. 

(ii) The boundary case. 

For the boundary Gaudin model ( [1V.5D we find that the eigenvalues are given by 



M 



E] = 2^{2zj) + aifizj + - E [V>izj - to) + ^{zj + t,) (IV. 10) 

with the constraints = 0, a = 1, 2, ■ ■ ■ , M, where 

t = c[ipit^-0 + ^{ta + 0]+2^{2to,) 

M N 

+ 2J2 - tp) + ^{ta + tp)] - ^[^{ta - Zk) + ^{ta + Z^)]. 

Pi^Q. k=l 

(IV.ll) 
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where a, c are some complex parameters. As previous section, the off-shell Bethe ansatz 
equations can be obtained by taking rj —> to the 



M 



a' 



(IV. 12) 



a=l 



Here t{zj),S{zj) and Fj^ can be found in $ , $^ are Bethe states given by 
= B{m° -2\ti)---B{m° -2M\zm)\0>, 

X B{m° - 2{a - l)|Wi) ■ ■ ■ B{m^ - 2M|tAf)|0 >, (IV. 13) 



where mP is a parameter and B is the (1,2) element of the 2x2 monodromy matrix 



which can be found in 23 



V Summary 

We have studied the elliptic Gaudin models and solutions of the corresponding KZ 
equations. In the first part of the paper, we have constructed the Gaudin models based 
on the face-type elliptic quantum groups. Hamiltonians of the models are diagonalized 
by using the algebraic Bethe ansatz. With the help of these Gaudin model Hamiltoni- 
ans, we have presented two types of elliptic KZ equations and give their solutions in 
the form of integrals. Then in the second part of the paper, we have constructed and 
diagonalized the boundary XYZ Gaudin model. 
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Note Added: After we finished this paper, we became aware that the integrable 
representation of solutions of the face-type elliptic KZ equations in the periodic case 

], ^ by means of techniques different from the ones 



have been obtained in 
presented in this paper. 



24 
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